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I. INTRODUCTION 



The bound-state-approach advocated by Callan and Klebanov(CK) ffl has been shown to 
work very well for static properties of strange baryons. In CK approach, strange baryons 
are described by properly quantized states of the i^-meson(s)-soliton bound system. Rho 
et al. H extended the CK approach further to baryons containing a heavy flavor such as 
charm or bottom. In particular, the mass spectra and magnetic moments M for charmed 
baryons are found to be strikingly close to the predictions of the quark model description. 
In these calculations, vector meson fields such as K*, D* and B* are eliminated in favor of 
a combination of a background and corresponding pseudoscalar fields, K, D and B. This 
approximation is valid only when vector mesons are sufficiently heavier than corresponding 
pseudoscalar mesons as in the case of p and ir {m p = 770 MeV, m n0 = 135 MeV: m^lm p = 
0.18). For charmed mesons or bottom flavored mesons, however, vector mesons are only 
a few percent heavier than corresponding pseudoscalar mesons: m K *, = 892 MeV, m K0 = 
498 MeV {m K0 /m K * = 0.56), m D , = 2010 MeV, m D0 = 1865 MeV {m D0 /m D , = 0.93), and 
m B * = 5325 MeV, m B0 = 5279 MeV (m B0 /m B , = 0.99). Thus it is needed to treat heavy 
vector mesons correctly on the same footing as heavy pseudoscalar mesons. 

The heavy quark symmetry is a new spin and flavor symmetry of QCD in the limit of 
infinite heavy quark masses. As a heavy quark becomes infinitely heavy, the dynamics of 
a heavy quark in QCD depends only on its velocity and is independent of its mass and 
spin. This symmetry can be seen in weak semileptonic decays |J, mass splittings and 
partial decay widths |J of heavy mesons and heavy baryons, whose masses are much bigger 
than the QCD scale, A QCD . Recently, effective heavy meson Lagrangians which have both 
chiral symmetry and heavy quark symmetry have been constructed by several authors J?],[8],|9[] . 
Also, a lot of works on heavy baryons as skyrmions a la Callan-Klebanov have been reported 
10| , p!T|JT2] , p!3[| . In a series of papers ||10|| , Jenkins et al. investigated the binding of a heavy 



meson with a soliton using such an effective Lagrangian. Nowak et al. [0 studied the heavy 
quark symmetry in heavy baryon mass spectra in connection with the Berry's phase. Gupta 
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et al. [0] discussed roles of light vector meson degrees of freedom such as u and p. In these 
works, however, only the leading order terms in the inverse of the heavy quark mass have 
been considered. Also, bound heavy mesons are assumed to sit at the center of the soliton 
with their wavefunctions taken as 5-functions. As a result, heavy mesons appear to be too 
deeply bound and Eg and Eg are degenerate in mass. 

In order to investigate more realistic cases with hyperfme splittings, one needs to include 
next to leading order terms in I/tuq. In Ref. ||10|| , mass corrections are roughly estimated 



by including mass differences between heavy pseudoscalar mesons and heavy vector mesons, 
while keeping the 5-function-like wavefunctions. Although it may work well for bottom fla- 
vored baryons, we may have some doubts on the validity of such 5-function-like wavefunctions 
for charmed baryons: the finite mass corrections need to be included in the wavefunctions 
of heavy mesons, leading to different radial functions, though sharply peaked at the center 
of the soliton. In this paper, we attempt to establish a "smooth" connection between the 
CK approach for light baryons and the heavy-meson-soliton bound state approach for heavy 
baryons by clarifying the above-mentioned problems. In our calculation, heavy pseudoscalar 
mesons and heavy vector mesons are treated on the same footing and the next to leading 
order terms in I /rug, are incorporated properly. 

In the following section, we introduce a simple Lagrangian which is relevant to our pur- 
pose. Then, a soliton-heavy-meson bound state is found in Sec. Ill by solving the equations 
of motion for the classical eigenmodes of heavy mesons moving in the soliton background. 
In Sec. IV, we discuss the mass formula for heavy baryons containing a heavy quark. We 
also discuss the heavy quark symmetry breaking by the Wess-Zumino term in the heavy 
baryon mass spectra. Section V contains summary and conclusion. 

II. MODEL LAGRANGIAN 

In order to avoid any unnecessary complications, we work with a simple Lagrangian for 
the interaction of light Goldstone bosons with heavy mesons, which has the SU (2) L x SU (2) R 
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chiral symmetry and the heavy quark symmetry in the heavy mass limit. One may obtain 
such a Lagrangian from the Skyrme model Lagrangian by trimming away all the higher 
derivative terms or from the heavy quark effective Lagrangian by including the next to 
leading order terms in 1/tjiq. 

Up to a single derivative on the Goldstone boson fields, the most general chirally invariant 
Lagrangian densityB can be written in a form of |J 



(1) 



+f Q + $*t^$) + y Q er»{*%A x *; + $?A X Q^), 

where $ and $* are the heavy pseudoscalar and the heavy vector meson doublet siwith 
masses M$ and M$*, respectively. For example, in the case of charmed mesons, we have 
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The Lagrangian density for the Goldstone boson fields is 
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U = ^ 2 = exp 
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and f n being the pion decay constant. The "Skyrme term" with a dimensionless parameter 
e is included to stabilize the soliton solution. Here f Q and g Q are the $$*7r and $*$*7r 



1 One may improve the model Lagrangian by including terms with more derivatives on the Gold- 



stone boson fields and incorporating vector mesons such as p and U) [13]. 

2 Here, we adopt a different convention for $ and $* than that of Ref. ||. Our 3>(<I>j!j) corresponds 
to their 
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coupling constants. The vector and axial vector potentials and are defined in terms 
of £ as 



(0c) 



and the covariant derivative and the field strength are 



Under SU(2) L x SU(2) R chiral transformations, the fields transform as 



£ £' = L£/it = fcfflt (U->U' = LUI&), 



(2) 



where L and R are global transformations in SU (2) L and SU (2) R respectively and h is 
a special unitary matrix depending on L, R and the Goldstone fields. Furthermore, the 
Lagrangian is invariant under the parity operation 

U(f,t) -> Pf/p- 1 = W(-f,t), 

$(f, t) -> V^V- 1 = -$(-f, t), (3) 
= -$;(-?,*). 

Here, we have used the fact that pions and heavy mesons (both pseudoscalar mesons and 
vector mesons) carry negative intrinsic parity. 

We have four parameters in the Lagrangian to be fixed; the pion decay constant f n , the 
Skyrme parameter e and the coupling constants f Q and g Q . The pion decay constant, /„-, 
and the Skyrme parameter, e, are fixed by fitting the nucleon and delta masses in the 577(2) 



sector |]T4] . As for the heavy meson coupling constant /_ and g Q , little has been known 



except the upper bound ||15fl . Thus, we use the heavy quark symmetry as a guide line. We 
use the empirical masses for the heavy meson masses, M$ and M$*. In the heavy mass 
limit, we have M$ ~ M$* with the mass difference being of order I/tuq at most and the 



two coupling constants f Q and g Q become related to each other by PJTB 



f Q = 2M^g Q (4) 

due to the heavy-quark spin symmetry. Furthermore, g Q approaches a universal constant g 
due to the heavy-quark flavor symmetry. 

The nonrelativistic quark model estimate of g (—0.75) || is consistent with the experi- 
mental value (|g| 2 ;$ 0.5) measured via the D* decay width and the D* + — > D + tt° and 
D* + — > D°ir + branching ratios jl7| . 

One may determine the two coupling constants f Q and g Q from a low energy chiral 



theory. In Ref. |18|], a Lagrangian for the interactions of K and K* mesons with pions is 
derived on the basis of 577(3) chiral symmetry along the hidden gauge symmetry scheme. 
Comparing it with our Lagrangian, we get f Q /2M K * = — ~ —0.71, which is very close 
to the nonrelativistic quark model prediction. Although the $*$*7r term proportional to g Q 
is missing in Ref. one can find such a term among the homogeneous solutions of the 



Wess-Zumino anomaly equation. (See Ref. [plj for further details.) Using the vector meson 



dominance hypothesis and the empirical value on the gx*^ coupling constant (~ 6), we 



obtain g Q ~ —0.7 from the chiral Lagrangian of Ref. |T6 



III. SOLITON-HEAVY MESON BOUND STATE 

The Lagrangian density C M supports a stable SU{2) soliton solution of "hedgehog" type 

Uo(r) — exp(ir ■ fF(r)), (5) 

with 

F(0) = 7i and F(r) ™ 0. (|a) 

The above solution carries a nontrivial winding number due to its nontrivial topological 
structure identified as the baryon number 

B= ^ / d ' r e^TriUld^Uld^Uld.U,) 

fro sin2jP m 



7r Jo r 2 



and a finite mass 

fl,-nn sin 2 F 1 sin 2 F sin 2 F 



M, 



sol 



f°° r> , K sin b . 1 sin ^ .sin b . 9 . , 

41 /„ r * y (F + 2 — » + 2?— <— + 2F » • <&=> 



with F' = ^f. 

ar 



Now, our problem is to find the eigenmodes of the heavy mesons moving in the static 
potentials provided by the B = 1 soliton configuration fl5|) sitting at the origin; viz., 

V» = (V°,V) = (0,iv(r)rxr), 
A<* = (A , A) = (0,|(ai(r)f + a 2 {r)ff ■ r)), 

with 

sin 2 (F/2) 



(6) 



v(r) 



> 

r 



a i( r ) = an d CL2{r) = b . (7) 

The equations of motion can be read off from the Lagrangian ([!]): 

(D^ + M|)$ = f Q A»%, (8) 

for the pseudoscalar meson field $ and 

D^ u + Ml*$* v = -f Q A v <$> + g Q e^ xp A x ^ p (9) 

for the vector meson fields $*. 

The conjugate momenta to the meson fields $ and $* are 

n« = = ($«°) t - 9 Q e^A v (10) 

respectively, and we get similar equations for IF and U*^ . Since I1q vanishes identically, the 
$q cannot be an independent dynamical variable. We eliminate the complementary <3>q field 
by using Eq. (||) 



$ *° = -^(AIT* + l 2 9 Q ^ k A k ^), (11) 
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which results in a set of coupled equations 

$* = _n*t _ 9q a x$* + -^3(3 ■ n*t) + -^-3(1- (3 x <£*)), 

n*t = 3 x (Z> x $*) + M 2 $* + / Q A$ - x n*t - g 2 Q A x (Ax $*) (12) 
2^ c 



1 x D {/3 ■ n*t + g Q A ■ (3 x $*)} . 



where = V — V. 

— * 

In order to express the equations of motion only in terms of $ and $*, we use the fact 
that the $q field is of order 1/ttiq at most; viz., 

~ -^r D & 1 ~ 0(1/M«.)- (13) 
Keeping this leading order term leads us to the equations of motion 

$* = -2g Q A x <S* - D x (D x $*) - Aff.$* 

-/ Q i$+^.n (M) 

Because of the spin-isospin mixing in the hedgehog configuration of the classical back- 
ground, the equations of motion (§) and ([14]) are invariant only under the rotation by the 
grand spin K defined by K = S+I+L with S(I) being the spin (isospin) of the heavy mesons 
and L the orbital angular momentum. Thus, eigenmodes are classified by the quantum num- 
bers k, rrik and P (the parity, P = (— with i being the orbital angular momentum) 
as 

k,m k ,P ^g-j 

<S*(f,t) = £ < mfc ,p(r,t)^ jmfc (f), 

k,m k ,P,k 

where 3^fc,p,m fc and D4,p,m fe are the generalized spherical spinor and vector harmonics, respec- 
tively, and k is an index to label the possible vector spherical harmonics with the same k, 
irik and P. To avoid cumbersome notation, we will suppress the trivial indices k, and P 
of the radial functions as ip(r,t) and ip* K (r,t). 
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From now on, we will restrict our consideration to k p = | + states, which are expected 
to have at least one bound state. Since pseudoscalar mesons do not carry spin, we have only 



one spherical spinor harmonics with k 1 



i + . 



yi 



i r 



t ■ r\±- 



(16) 



Here, x± is the isospin basis for the heavy meson doublets, i.e. 











x+ = 




and X- = 













For vector mesons with spin 1, we can construct two different k p 



i + 



(17) 



vector spherical 



harmonics 19 



: viz., 



Putting 



2 '^'^ 2 



1 



rx±, 



'An 

y?\ ± i(f) =%— ={t x r) X ±. 

2'" t "' =c 2 V07T 



(18) 



- + ±- 

2 ' ' 2 



+^(r)e-*-*3^° +±1 (r) ) 



(19) 



(20) 



into the equations of motion (||) and ([14]), we obtain three coupled differential equations for 
the radial functions: 



<p" + -if' + (lu 2 - M, 



fc 



1 



(21) 



-)<p = 2v{v - ^)cp + -l{a x + a 2 )^i - -^f^tp*, 

<P? + %X + - Ml. - ^)<p\ = f -f{a x + a 2 )y? + 2v\\ 

+V2(g Q a 1 u --v + v')<p* 2 , 

2 2 / r 1 

^2" + -V2 + _ M l* 2 ")v 9 2 = — % a i¥ + v2(a;y Q ai v + v')tp\ 

r r z a/2 r 

4 

+ (-^ Q (ai + a 2 ) v + 4w 2 )^. 

The wavefunctions are normalized such that each mode carries one corresponding heavy 
flavor number: 



1 = ^dr {2lu [\ V \ 2 + |y?*| 2 + |^| 2 ] + 9 Q [(a, + a 2 )\<p* 2 \ 2 - y/2 ai (<p? <p* 2 + ^V*i)] } , (22) 

where we have kept terms up to the next to leading order in 1/toq. 
Near the origin, the equations of motion behave asymptotically as 







2 . 

+ V = o, 

r 




H>\ + 


r 


V- 


2^ , 

r 2 ^2 


^2 + 


2 <n*> 


2 * 

^2 = 


2v/2 , 

r 2 



(23) 



They imply that we have three independent solution sets as 

(a) v?(r)=v?(0)+O(r 2 ), 
tf(r) = 0(r 2 ), 

(b) v?(r)=0(r 2 ), 

^*(r) = <(0) + O(r 2 ), 

(c) v?(r)=0(r 4 ), 

^W = ^'(0)r 2 + O(r 4 ), 



(24) 



with V2ip* bl (0) = </4(0) and ^*i(0) = -y/2(p%(0). For sufficiently large r(> 1/M*), the 
three equations decouple from each other: for example, 

2 



^" + V + (w 2 -M!»)y? = o. 

Thus the bound state solutions (uo < M$) are 

</?(r) = a- 



(25) 



(pl(r) = ai- 



^j(r) = a 2 - 



r 



(26) 



with three constants a, cti and a 2 - 
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The lowest energy bound states are found numerically, and the results are shown in 
Table I and Fig. 1. In Table I, the input parameters are listed together with the numerical 
results on the lowest bound states. In Fig. 1, we give the radial functions <£>(r) and <fl(r) 
for the D and D* mesons (solid curve) and the B and B* mesons (dashed curves). By 
comparing the two cases, one can easily check that as the meson mass becomes larger, (1) 
the radial function becomes more sharply peaked at the origin and (2) the role of the vector 
mesons becomes important so that the radial function (fl(r) becomes comparable to <£>(r) 
(see also the ratio (pl(0)/(p(0)). The radial function ^(r), though not shown in Fig. 1, is 
hardly distinguishable from \f2Lp\(r). This can be understood as follows: due to their heavy 
masses, heavy mesons are localized in the region r^l/M$, where 

Mr) + a a (r)] ~ [-^(r)] ~ F'(0) + 0(r 2 ), 

1 (27) 
t,(r)~--±F /2 (0)r + ---, 

so that the equation of motion for (ip\ — ^y 9 ^) is completely decoupled from those for ip and 

It would be interesting to compare our radial functions with those of Ref. || and Ref. 
||1 0|| . In Ref. P|, vector mesons are assumed to be sufficiently heavy and the following ansatz 
is made: 

•S = («) 

which implies that 

¥>l(r) = ^ M (ai(r) + a 2 {r))<f(r), 

ti(r) = --^- ai (7>(r). (29) 

As v^Vi ~ ^2 f° r heavy mesons due to Eq. fl2T|), we have only to compare Lp\ with ip 
in Eq. (p5|). In the heavy mass limit, both should play equally important roles. But the 
ansatz strongly suppresses the role of vector mesons by a factor of y/^ef^/M^*, since one 
obtains F f (0) ~ —2ef n in the Skyrme-term-stabilized soliton solution. For example, this 
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factor amounts to 0.56, 0.25 and 0.09 for the cases of M K *(892 MeV), M D *(2010 MeV) and 
Ms* (5325 MeV), respectively. Therefore, the ansatz of Eq. ( p8| ) is not valid unless the 
vector meson is much heavier than the corresponding pseudoscalar meson. 



The wavefunctions of Refs. [10,T6fl are obtained in the heavy mass limit, M $ , M^, — > oo 



and can be written in our convention as 



l l 



(30) 



where the radial function f(r), normalized as Jr 2 dr\f\ 2 = 1, is strongly peaked at the origin. 
It implies that 

¥>(r) = -^(r) = -i?S(r) ~ 1^=/^)- W 



These radial functions satisfy the normalization condition of Eq. (22) in the leading order 
in l/m Q ; viz., 

POO 

2lu b r 2 dr(|v9| 2 + |^| 2 + |^| 2 ) = l. (0b) 

J 

It is interesting to note that the pseudoscalar meson and three vector mesons contribute 
equally to the bound state. 

Comparing our numerical results given in Table I with the binding energy Eb = 



— ~g Q F'(0) of Refs. fl0l , [T6"|l which gives ~ 800 MeV with the same input parameters, one can 
see that the 1 / tuq corrections amount to ~ 200 MeV in the bottom sector and ~ 300 MeV 
in the charm sector. This is one of the main results of this work. 



In Ref. |I(J, the rms radii of the heavy flavor current in heavy baryons are essentially 
zero. Due to the 1/mg corrections, however, we have non-zero finite size rms radii in our 
calculation, viz. ~ 0.3 fm for bottom flavored baryons and ~ 0.4 fm for charmed baryons. 
This implies that the rms radii of heavy flavored baryons become small as the masses become 
large. Due to this effect, the binding energy is smaller than the one obtained with 5-function- 
type solutions. 
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IV. HEAVY BARYONS AND HYPERFINE SPLITTINGS 



So far we have considered soliton-heavy-meson bound states to the order N® with N c 
being the number of color. The combined system of the soliton and a bound heavy meson 
carries a baryon number and a heavy flavor number, but does not have the spin and isospin of 
a heavy baryon. Up to order N®, the soliton-heavy-meson bound state should be understood 
as a mixed state of three degenerate heavy baryons containing a heavy quark Q; Eq, Aq 
and Eq, whose mass is M so i + ujb- In order to give the spin and isospin quantum numbers 
and the hyperfme splittings, we have to go to the next order in 1/N C , i.e., C^A^ 1 ). This is 
done by quantizing the zero modes associated with the simultaneous SU(2) rotation of the 
combined system. A standard collective coordinate quantization procedure leads us to the 
mass formula for a heavy baryon with spin J and isospin /: 

M = M sol + uj b 

+^ (cJ(J + 1) + (1 - c)I(I + 1) + - A c{c - 1)) 

+0(1/M 2 ). (31) 

Here X is the moment of inertia of the soliton configuration against the SU(2) collective 
rotation: 

and c is the hyperfme splitting constant which can be obtained by directly applying the 
techniques developed in Ref. [pTSfl : 



roo 

c = r 2 dr { 2uji 



o 

+ k 



(M 2 - l M? - ^2l 2 ) - f cos 2 (F/2)(|yf - |^| 2 ) 



{F > _ S ^l)\^ - i=^(3cosF + 1)(^V 2 + rffc) 
r y/2 r 



(0b) 



3^Q 

We note that we have also kept terms of the next to leading order in 1/ttiq. One can easily 
see that the hyperfme constant c is of order l/rnq. The leading order terms proportional to 
uj b vanish identically when the radial functions of Eq. (|30|a) are used. 

According to the formula, the masses of heavy baryons containing a single heavy quark 
have following hyperfme splittings: 
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By eliminating c from Eq. (|32|) we have a model independent relation 

\{2M Vq + M SQ ) - M Aq = |(M A - Mjv). (33) 

With the experimental values M^ p '(= 2453 MeV), M^ p '(= 2285 MeV) and M^ p '(= 
5641 MeV), we predict the mass of S* to be 2493 MeV and the averaged mass M^ b (= 
|(2Me* + MeJ 5836 MeV. Since c is of order 1/rriQ, the masses of £q and of Eg are 



degenerate in the infinite mass limit as the heavy quark symmetry implies and Eq. (|33"D is 
reduced to M Sq - M Aq = |(M A - M N ) as in Refs. ■ 

Numerical results (Result I) on the heavy baryon masses are shown in Table II. They 
are in rough agreement with the experimental values. Result II is obtained by taking the 
two coupling constants as free parameters. To fit the experimental masses of A c and E c , 
one should have f Q /2Mo* = —1.04 and g Q = —0.40, which implies that the heavy quark 
symmetric relation is strongly broken in the charm sector. Note that as far as the two 
coupling constants are related by Eq. (||), the hyperfine constant is too small. 

In order to improve the situation, one may consider higher order terms in the 1 /rrtq ex- 
pansion or higher derivative terms of the pion fields. As a guide line, we may use the Skyrme 



Lagrangian |L8| with the vector mesons included via the hidden gauge symmetry, since in 
the strangeness sector the heavy quark symmetry becomes no longer a good symmetry but 
the 577(3) chiral symmetry becomes rather a good symmetry. 

Among many possible terms which will be discussed below, the Wess-Zumino (WZ) term 
is known to play the most important role in CK approach PJ: 

Avz = -^p' i&D^ - (D^*) . (34) 

Here, N c is the number of color, is the $-meson decay constant and is the topological 
baryon number current. Although its role fades out in the heavy mass limit, the WZ term 
should not be disregarded in the case of finite quark masses. 
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In addition, there are other contributions of order I/tuq to the binding potential. We 
introduce a typical A ■ A potential in the Lagrangian as in CK approach, which is next to 
leading order in the derivative of pion fields and turns out to have non-negligible effects in 
the strangeness sector: 

£ (2) = -&A^$. (35) 
Now we discuss the effects of the above terms in detail. Let us write 

5£ = £ wz + £(2)- (36) 

This additional Lagrangian modifies the equations of motion for the pseudoscalar meson 
field $ as 



(D^ + Ml)* = f Q A^l - —p^B.D^ 



2iN c 
4/1 

(37) 



while those for the vector meson field $* remain the same as Eq. ([14]). Consequently, the 
radial function of the k p = ^ + eigenmodes is altered as 



ip" + V + (lu 2 -Ml- hip 

2 

= |/ Q (ai + a 2 )0i - ^/qCIi^ - (2^A - 2v{v ) - \(3aj + 2a 1 a 2 + a\))ip 



where 

N c 1 sin 2 F . 

Hr) = - Ti ^— (39) 

The WZ term contributes to the hyperfine splitting constant c as 

/•OO 

Sc = 2 / r 2 dr\p\ 2 \, (40) 



and to the normalization condition of Eq. (|22| ) by the same amount. Note that there is no 
direct contribution from £( 2 ) to this quantity. 

We begin with the role of 5£ in the strangeness sector, where the above model lagrangian 
does not work well as expected in this sector. In Table III, we show the numerical results for 
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strange baryons obtained with the input parameters f n = 64.5 MeV, e = 5.45, Mr- = 495 
MeV, M K , = 892 MeV and f Q /2M K * = -4- = g Q . Note that the role of £ (2) in the binding 
energy is important (~ 80 MeV), while its effect on the hyperfine constant c is rather small. 
The Wess-Zumino term plays a crucial role in the hyperfine constant c, of which more than 
80% comes from the WZ term. As shown in Fig. 2, the effect of the WZ term on the radial 
wavefunction is also remarkable; with the WZ term, the vector meson contribution to the 
bound states is much suppressed compared with that of the pseudoscalar meson. 

In the charm sector the role of the Wess-Zumino term in the heavy mass limit is weakened 
as discussed in Ref. ||11|| . This results from the fact that the role of vector mesons balances off 



that of pseudoscalar mesons in the heavy mass limit. The decoupling of the WZ term in the 
heavy mass limit is originally argued in Ref. pl~| , p2| . In the chiral limit where pseudoscalar 
mesons predominate, the WZ term is entirely expressed in terms of these pseudoscalar 
mesons. To take into account the WZ term in the finite mass region, we consider the most 
characteristic expression of the WZ term p|fl8[| constructed in terms of pseudoscalar mesons 
and an adjustable parameter 7. The parameter 7 contains the trace of cancellation between 
contributions of the vector and of the pseudoscalar mesons and should depend on 1/rriQ. 
That is, we take 

5C = 1 C wz + eC {2) (41) 

with the same £ wz and Cr 2 ) as given by Eqs. (|34|,|35]). The parameter e has the role of 
turning on and off the effect of £( 2 ). Here, /$ is the D-meson decay constant f D , which is 
known to be 1.8 times larger than the pion decay constant Although the C( 2 ) plays a 
minor role for the heavy flavors such as charm, we keep it to compare its effects in the charm 
sector with those in the strangeness sector. In Fig. 3, we present uj b and c as a function of 
the mitigating factor 7. The role of £( 2 ) is shown as narrow stripes, with ~ 30 MeV effect 
on the energy and ~ 0.04 on the hyperfine constant. However, as we can see in Fig. 3, the 
dependence of the mass spectrum on the parameter 7 is not negligible. In order to fit the 
charmed baryon masses, we need to have uj b = 1416 MeV and c = 0.16. Then we have 

16 



M Ac = 2285 MeV 

M Sc = 2449 MeV, (42) 
M s * = 2495 MeV. 

Also from Fig. 3, one can estimate the mitigating factors as 7 ~ 0.25 and e = 1, which 
reproduce the above mass spectra and then we have y (r 2 ) c = 0.37 fm. It implies that the 
increase in /$ alone is not enough to take fully into account the role of the Wess-Zumino 
term in the charm sector. 

The dependence of 7 on meson masses can be derived by showing how the Wess-Zumino 
term scales out as the mass increases. We are not in a position to illustrate this dependence 
yet. However, if we assume the dependence to be inversely proportional to the m$, we find 
7 ~ 0.5 GeV/m$. Note the coincidence of the 7-factor (~ 0.25) with the meson mass ratio, 
m K /m D . 



V. SUMMARY AND CONCLUSION 

In this work, we have investigated the mass spectrum of heavy baryons containing a 
single heavy quark in the bound state approach of Skyrme model. To this end, we have 
worked with the heavy meson Lagrangian of Ref. || which includes the I/tuq order terms. 
The large binding energy obtained in the infinite mass limit is lowered by introducing 1 /tuq 
corrections. The binding energy is changed from ~ 800 MeV to ~ 500 MeV for D(D*) 
mesons and to ~ 600 MeV for the B(B*). The effect may be crucial for the loosely bound 
exotic states such as "pentaquark" baryons [^3[^4|. However, due to the realization of the 



heavy quark symmetry, the hyperfine splitting constant comes out too small compared with 
the experimental one. For example, we get c = 0.05 for charmed baryons while it should 
be ~ 0.14 to reproduce the experimental masses. To resolve this problem, we introduce 
the WZ term in a mitigated form, which is known to have a crucial role in the strangeness 
sector. To reproduce the experimental masses for charmed baryons, its strength should be 
weakened by a factor of 4. 
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FIGURES 

FIG. 1. ip(r) and ip*(r) for Q = c (solid) and b (dashed), <p%{ r ) is nearly equal to v / 2</ ? i( r ) f° r 
both cases. 

FIG. 2. ip(r) and <p\{r) for (a) B and B*, (b) D and D*, (c) K and K* with (solid) and without 
(dashed) the Wess-Zumino term. Each fields are normalized as / drr 2 \ip\ 2 = 1. 

FIG. 3. oj b and c vs. 7 obtained for charmed baryons. 
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TABLES 



TABLE I. Summary on the input parameters and the numerical results on the bound state 



Q 


,a) 
J-rr 


e c ) 


M$ 




fa) 
J Q 




a) 
^B 




c c ) 


tf(o)Mo) 


c 

b 


64.5 
64.5 


5.45 
5.45 


1872 
5275 


2010 
5325 


-3016 
-7988 


-0.75 
-0.75 


1481 
4722 


0.39 
0.29 


0.05 
0.02 


-0.828 
-0.932 








a ) in MeV unit, ^ in fm unit, and c ) 


dimensionless quantities. 










TABLE II. 


. Numerical result 


s on the heavy baryon masses. 




Q 






/ Q /2M*. 




9q 


a) 


c 


A Q 




a 


c 


exp. 
I 

II 




-0.75 
-1.04 




-0.75 
-0.40 


1481 
1419 


0.05 
0.14 


2285 
2348 
2287 


2453 
2535 
2454 


2548 
2497 


b 


exp. 
I 


6 ) 


-0.75 




-0.75 


4722 


0.02 


5641 
5589 


5781 


— 
5786 










a) i n 


MeV unit, b ) Particle Data Group @. 












TABLE III. 


WZ term and Strang 


;e baryon masses. 










£ (2) 


a) 




c 


M A 




M$ 


M& 




off 
off 
on 
on 


exp. 


off 
on 
off 
on 


389 
291 
191 
109 




0.098 
0.148 
0.717 
0.791 


1257 
1160 
1095 
1022 
1116 




1433 
1326 
1151 
1063 
1192 


1462 
1369 
1361 
1295 
1385 


0.62 
0.56 
0.41 
0.39 



) in MeV unit, ^ in fm unit. 
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Fig.2 : tp(r) and ip\(r) for (a) B and B*, (b) D and D*, (c) K and K* 
with (solid) and without (dashed) the Wess-Zumino term. 
Each fields are normalized as / drr 2 \(p\ 2 = 1. 
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Fig. 3 : uj b and c vs. 7 obtained for charmed baryons. 
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